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1. Introduction 

In deriving uniform asymptotic expansions of a certain class of integrals one encounters 
the problem of expanding a function, that is analytic in some domain f2 of the complex 
plane, in two points. The first mention of the use of such expansions in asymptotics is 
given in [1], where Airy-type expansions are derived for integrals having two nearby (or 
coalescing) saddle points. This reference does not give further details about two-point 
Taylor expansions, because the coefficients in the Airy-type asymptotic expansion are 
derived in a different way. 

To demonstrate the application in asymptotics we consider the integral 

J?iH = i/ c eU( ^ )/(z)(l2 ' (1) 

where w is a large positive parameter and b is a parameter that may assume small 
values. The contour starts at ooe~ l ' K / 3 and terminates at ooe l7r / 3 , and lies in a domain 
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where the function / is analytic. In particular, / is analytic in a domain that contains 
the saddle points ±6 of the exponent in the integrand. One method for obtaining an 
asymptotic expansion of Fb(uj) that holds uniformly for small values of b is based on 
expanding / at the two saddle points: 

oo oo 

Hz) = Mz 2 - b 2 r +zY,B n (z 2 - b 2 r, (2) 

n=0 n=0 

and substitute this expansion into (1). When interchanging summation and integration, 
the result is a formal expansion in two series in terms of functions related with Airy 
functions. A Maple algorithm for obtaining the coefficients A n and B n , with applications 
to Airy-type expansions of parabolic cylinder functions, is given in [4] . 

In a future paper we use expansions like (2) in order to derive convergent expansions 
for orthogonal polynomials and hyper geometric functions that also have an asymptotic 
nature. The purpose of the present paper is to give details on the two-point Taylor 
expansion (2), in particular on the region of convergence and on representations in 
terms of Cauchy-type integrals of coefficients and remainders of these expansions. Some 
information on this type of expansions is also given in [6], p. 149, Exercise 24. 

Without referring to applications in asymptotic analysis we include analogous prop- 
erties of two-point Laurent expansions and of another related type, the two-point Taylor- 
Laurent expansion. 



2. Two-point Taylor expansions 

We consider the expansion (2) in a more symmetric form and give information on the 
coefficients and the remainder in the expansion. 

Theorem 1. Let f(z) be an analytic function on an open set Q C (D and z\, z 2 € 
with Z\ 7^ z 2 . Then, f(z) admits the two-point Taylor expansion 

N-l 

f( z ) = i a n(zi,z 2 )(z - zi) + a n (z 2 , zi)(z - z 2 )\ (z-z 1 ) n (z-z 2 ) n +r N (z 1 ,z 2 ;z), (3) 

n=0 

where the coefficients a n (zi,z 2 ) and a n (z 2 ,zi) of the expansion are given by the Cauchy 
integral 

( \ — * f f( w ') ^ w (A\ 

MZl ' Z2) = 2^{z 2 - zV) Jc (w- Zl nw-z 2 )^- W 

The remainder term r^{z\, z 2 \ z) is given by the Cauchy integral 

1 I' f(w)dw _ N N 



r N ( Zl ,z 2 ;z) = — — p- — Az-zi) (z - z 2 ) . (5) 

2m J c (w — zi) n (w — z 2 ) N (w — z) 

The contour of integration C is a simple closed loop which encircles the points z\ and 
z 2 (for a n ) and z, Z\ and z 2 (for r N ) in the counterclockwise direction and is contained 
in f2 (see Figure 1 (a)). 
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The expansion (3) is convergent for z inside the Cassini oval (see Figure 2) 
ZuZ2 = {zen, \{z - Zl )(z - z 2 )\ < r} 

where 

r = lnf we€ \ n {\(w - zi)(w - z 2 )\} . 

In particular, if f(z) is an entire function (fi = (D), then the expansion (3) converges V 
z £ <D. 

Proof. By Cauchy's theorem, 



1 f f(w) dw 



lc w - z 

where C is the contour defined above (Figure 1 (a)). Now we write 

1 Z + W — Z\ — Z 2 1 



(6) 



where 



w — z (w — zi)(w — z 2 ) 1 — u ' 
(z - zi)(z - z 2 ) 



u = 



(w - z 1 )(w - z 2 ) ' 



(7) 



(8) 



Now we introduce the expansion 



^ l-u 



1-u 



(9) 



n=0 



in (7) and this in (6). After straightforward calculations we obtain (3)- (5). 

For any z £ Zl>Z2 , we can take a contour C in O such that \(z — z\)(z — z 2 )\ < 
\(w — z\){w — z 2 )\ V w £ C (see Figure 1 (b)). In this contour is bounded 

by some constant C: \f(w)\ < C. Introducing these two bounds in (5) we see that 
limAT^oo rjv(zi, z 2 ; z) = and the proof follows. □ 





(a) (b) 
Figure 1. (a) Contour C in the integrals (3)-(5). (b) For z G O 

zi,z2> we can take a contour 

C in £1 which contains ZljZ2 inside and therefore, \(z — Zi)(z — z 2 )\ < \(w — Z\)(w — z 2 )\ V 
w e C. 
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Figure 2. Shape of the Cassini oval depending on the relative size of the parameter r and the 
focal distance \zi — z 2 \ ■ 



2.1. An alternative form of the expansion 

The present expansion of f{z) in the form (3) stresses the symmetry of the expansion 
with respect to z\ and z 2 . In this representation it is not possible, however, to let 
z\ and Z2 coincide, which causes a little inconvenience (the coefficients a n (zi,z 2 ) be- 
come infinitely large as z\ — > z 2 ; the remainder rj^{z\^z 2 \z) remains well-defined). An 
alternative way is the representation (cf. (2)), 

00 

f(z) = [M*i, **) + B n ( Zl ,z 2 ) z) (z - z 1 ) n (z - z 2 ) n , 

n=0 

and we have the relations 

A n (z 1 ,z 2 ) = -z 1 a n (z 1 ,z 2 ) - z 2 a n (z 2 ,z 1 ), 
B n (z 1 ,z 2 ) = a n (z 1 ,z 2 ) + a n (z 2 , zi), 

which are regular when z\ — ► z 2 . In fact we have 

A I \ 1 f w- z x - z 2 
A n (z 1 ,z 2 ) = — / - — T f{w)dw, 

r ( \ - — f /(^) dw 

^{ZuZ2)- 27r .J c[{w _ Zi){w _ Z2)]n+1 . 

Letting z\ — > and z 2 — ► 0, we obtain the standard Maclaurin series of f(z) with even 
part (the A n series) and odd part (the B n series). 

2.2. Explicit forms of the coefficients 

Definition (4) is not appropriate for numerical computations. A more practical formula 
to compute the coefficients of the above two-point Taylor expansion is given in the 
following proposition. 

Proposition 1. Coefficients a n (zi,z 2 ) in the expansion (3) are also given by the for- 
mulas: 

a (z 1 ,z 2 ) = ^- (10) 

z 2 — Z\ 
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and, for n = 1, 2, 3, 



_ " („ + fe _ i)i (-l)"+i n /("- fc )(z 2 ) + (-l) fc fc/("- fc )( Zl ) 



n\(zi - z 2 ) 



Proof. We deform the contour of integration C in equation (4) to any contour of the 
form C\ U C 2 also contained in U, where C\ (C 2 ) is a simple closed loop which encircles 
the point z\ (z 2 ) in the counterclockwise direction and does not contain the point z 2 (z\) 
inside (see Figure 3 (a)). Then, 



a n (z 1 ,z 2 ) 



2m(z 2 - zi) \J Cl (w - z 2 ) n+1 (w 
1 f 1 d n ~ x f(w 



dw f 



dw 



(z 2 ~ zi) I (n - 1)! dw n ~ l (w - z 2 ) n+1 



c 2 {w - Zl ) n {w - z 2 ) n+1 
1 d n f{w) 



+ 



n\ dw n (w — zi) T 



From here, equations (10)-(11) follows after straightforward computations. 



□ 






(a) (b) (c) 

Figure 3. (a) The function (w — 2 2 ) _n_1 /(w) is analytic inside C\, whereas (w — Zi)~ n f(w) 
is analytic inside C 2 . (b) The function (w — z 2 )~ n g\ (w) is analytic inside C\, whereas (w — 
Z\)~ n g 2 (w) is analytic inside C 2 . (c) The function (w — z 2 )~ n ~ 1 g(w) is analytic inside C\, 
whereas (w — Zi)~ n f(w) is analytic inside C 2 . 



2.3. Two-point Taylor polynomials 

Next we can define the two-point Taylor polynomial of the function f(z) at in the 
following way: 

Definition 1. Let z be a real or complex variable and z\ and z 2 (z\ 7^ z 2 ) two real or 
complex numbers. If f{z) is n — 1— times differ entiable at those two points, we define 
the two-point Taylor polynomial of f(z) at z\ and z 2 and degree 2n — 1 as 

n-l 

P n (z 1 ,z 2 ;z) = ^2 [ a k{zi,z 2 )(z - zi) + a k (z 2 ,z 1 )(z - z 2 )] (z - zi) k (z - z 2 ) k , 

k=0 

where the coefficients ak(zi,z 2 ) are given in (10)-(11). 
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Proposition 2. In the conditions of the above definition, define the remainder of the 
approximation of f(z) by P n (z\, z 2 ; z) at z\ and z 2 as 

r n (z 1 ,z 2 ;z) = f(z) - P n {z 1 ,z 2 ;z). 

Then, (i) r n (z\, z 2 ;z) = o(z — z\) n ~ l as z — ► z\ and r n (z\, z 2 ; z) = o(z — z 2 ) n ~ l as z — > 
z 2 . (ii) If f(z) is n— times differentiable at z\ and/or z 2 , then r n (z\, z 2 ; z) = 0{z — z\) n 
as z — ► z\ and/or r n (z\, z 2 ; z) = 0{z — z 2 ) n as z — > z 2 . 

Proof. The proof is trivial if f(z) is analytic at z\ and z 2 by using (5). In any case, for 
real or complex variable, the proof follows after straightforward computations by using 
l'Hopital's rule and (10)-(11). □ 
Remark 1. Observe that the Taylor polynomial of f(z) at z\ and z 2 and degree 2n — 1 
is the same as the Hermite's interpolation polynomial of f(z) at z\ and z 2 with data 
f(z l ),fi(z l ),...,f( n - 1 \z l ),i = 1,2. 



3. Two-point Laurent expansions 

In the standard theory for Taylor and Laurent expansions much analogy exists between 
the two types of expansions. For two-point expansions, we have a similar agreement in 
the representations of coefficients and remainders. 

Theorem 2. Let Q and $7 be closed and open sets, respectively, of the complex plane, 
and r2 C SI C (D. Let f(z) be an analytic function on SI \ Sl an d z\, z 2 G S7 with 
Z\ z 2 . Then, for any z6fi \ S7 ; f(z) admits the two-point Laurent expansion 

N-l 

f(z) = K(zi,z 2 )(z - zi) + b n (z 2 , Zi)(z - z 2 )\ (z - Zl ) n {z-z 2 ) n + 

n=0 

N-l ( 12 ) 

S ^[c n {z l ,z 2 ){z-z l )+c n {z 2 ,z l )(z-z 2 )](z-z l ) 11 1 (z-z 2 ) n 1 + 

n=0 

rN{zi,z 2 ;z), 

where the coefficients b n (zi,z 2 ), b n (z 2 ,z\), c n (z\,z 2 ) and c n {z 2 ,z\) of the expansion are 
given, respectively, by the Cauchy integrals 

b n {z 1 ,z 2 ) = —- -/ — — — 13 

2m(z 2 - zi) J r (w - zi) n {w - z 2 ) n+i 



and 



C n (zl,Z 2 ) = 1 f (W - Zl ) n+1 { W - Z 2 ) n f(w)dw. 

2m(z 2 - z\) J F2 

The remainder term r^(zi, z 2 ; z) is given by the Cauchy integrals 

( \ _ 1 [ f{w)dw N N 

r N ( Zl ,z 2 ;z) =— Az-zi) \z-z 2 ) - 

2m J Vi (w — zi) n (w — z 2 ) N (w — z) 



(14) 



-/ 

2m J r 



(w — Z\) (w — z 2 ) f(w)dw 1 

v7^~Z ( Z - Zl )N( z _ Z2 )N- 



(15) 
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In these integrals, the contours of integration T\ and F 2 are simple closed loops contained 
in which encircle the points z\ and z 2 in the counterclockwise direction. Moreover, 

F 2 does not contain the point z inside, whereas Fi encircles F 2 and the point z (see 
Figure 4 (a)). 

The expansion (12) is convergent for z inside the Cassini annulus (see Figure 5) 
A ZuZ2 = {zen\n , r 2 < \{z - Zl )(z - z 2 )\ <n} (16) 

where 

n = Inf weC \ n {\(w - zi)(w - z 2 )\] , r 2 = Sup^ Gf7n {\{w - zi)(w - z 2 )\} . 



Proof. By Cauchy's theorem, 

1 J f(w)dw 1 I f(w) dw 

where Ti and F 2 are the contours defined above. We substitute (7)-(8) into the first 
integral above and 

1 Zi+Z 2 — z — w 1 _ (w — Zi)(w — z 2 ) 



w — z (z — zi)(z — z 2 ) 1 — u (z - zi)(z — z 2 ) 

into the second one. Now we introduce the expansion (9) of the factor (1 — u) _1 in both 
integrals in (17). After straightforward calculations we obtain (12)-(15). 

For any z verifying (16), we can take simple closed loops Ti and F 2 in Q\Q such that 
\{z — z\)(z — z 2 )\ < \(w — zi)(w — z 2 )\ y w G Ti and \ (z — z\)(z — z 2 )\ > \{w — z\)(w — z 2 )\ 
V w G F 2 (see Figure 4 (b)). On these contours |/(w)| is bounded by some constant 
C: \f(w)\ < C. Introducing these bounds in (15) we see that liniAr^oo tn{z\, z 2 ; z) = 
and the proof follows. □ 





(a) (b) 
Figure 4. (a) Contours Fi and F 2 in the integrals (12)-(15). (b) For z G A Zl:Z2 , we can take 
a contour F 2 in Q situated between Qq ar, d A ZljZ2 and a contour F\ in £1 which contains A ZlyZ2 
inside. Therefore, \(z — Z\){z — z 2 )\ < \(w — Zi){w — z 2 )\ V w £ Fi and \(w — Z\){w — z 2 )\ < 
\{z - Zl )(z - z 2 )\ VweF 2 . 



(a) 4n > 4r 2 > \z\ - z 2 \ 2 (b) 4n > \z\ - z 2 \ 2 = 4r 2 (c) 4n > |zi - z 2 | 2 > 4r 2 





(d) 4ri = |zi — z 2 | 2 > 4r 2 (e) \zi — z 2 \ 2 > 4j"i > 4r 2 

Figure 5. Shape of the Cassini annulus depending on the relative size of the parameters r\, r 2 
and the focal distance \z\ — z 2 \. 

If the only singularities of f(z) inside S7q are just poles at Z\ and/or z 2 , then an 
alternative formula to (13) and (14) to compute the coefficients of the above two-point 
Laurent expansion is given in the following proposition. 

Proposition 3. Suppose that gi(z) = (z — zi) mi f(z) and g 2 (z) = (z — z 2 ) m2 f(z) are 
analytic functions in Q for certain m\, m 2 G IN. Then, for n = 0,1,2,..., coefficients 
b n (zi,z 2 ) and c n (z\,z 2 ) in the expansion (12) are also given by the formulas: 



. , n _ fn + mi-n (-l) fc+1 (n + l) k g\ ' ( Zl ) 

AU 2 '~ ^ \ k ) ( n + mi -l)\( Zl - Z2 )n+k+2 + 

' ~- + m 2 \ (-l) k (n) k g 2 + \z 2 ) 



k = 

n 



E 



k=0 



(18) 



k J (n + m 2 )\(z 2 - z 1 ) n+k + 1 ' 
where (n)k denotes the Pochhammer symbol and 

m 1 —n — 2 , v / \ / \„_J._] (mi— n — k — 2)/ \ 

'mi - n - 2 \ / n \ (z x - z 2 ) n k l g\ ' (z x ) 



c n (z u z 2 )=- J2 kl ( mi k 2 )( n k ) 

k=o v v ■' 

-2—n — l 



(mi — n — 2)! 



/ m 2 - n - 1 \ / n+ 1 \ (z 2 - z x ) n k g y 2 2 \z 2) 



(19) 



fc=0 



J \ k J (m 2 — n—l)\ 



In these formulas, empty sums must be understood as zero. Coefficients b n (z 2 ,Z\) and 
c n (z 2 ,Zi) are given, respectively, by (18) and (19) interchanging Z \, g\ and mi by z 2 , 
g 2 and m 2 respectively. 
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Proof. We deform both, the contour T 1 in equation (13) and T 2 in equation (14), to 
any contour of the form C\ U C 2 contained in $7, where C\ (C 2 ) is a simple closed loop 
which encircles the point z\ {z 2 ) in the counterclockwise direction and does not contain 
the point z 2 (zi) inside (see Figure 3 (b)). Then, 



b n (z 1 ,z 2 ) 



1 



2ni(z 2 - zi) 
92(w) 



9i(w) 



dw 



dw 



C2 (w - zi) n (w - z 2 ) n + m2+1 



jn+mi-l 



n+1 (w - z x ) n + r - 



+ 



z i — z i I ( n + fn\ — 1)! dw n ~ 
1 d n+m2 g 2 (w) 



(w - z 2 ) n+1 



+ 



W = Z\ 




(n + m 2 )\ dw n+ni2 (w-ZxY 
(w - z 2 ) n g 1 (w" 



2m(z 2 -z 1 ) [J Cl (w - zi)" 1 !""" 1 

dmi -n-2 \( W - Z2 y gi ( yW 



dw + 



dwmi -n-2 y ( mi _ n _ 2)1 



+ 



IC 2 
d m 2 — n— 1 



(w - zi) n+1 g 2 (w) 
(w - z 2 ) m2 ~ n 



-n-1 



(W ~ -2l) n+1 ff 2 (w) 

(m 2 — n — 1)! 



□ 



From here, equations (18) and (19) follow after straightforward computations 

Remark 2. Let z be a real or complex variable and z\,z 2 (z\ / z 2 ) two real or 
complex numbers. Suppose that gi(z) = (z — zi) mi f(z) is n— times differentiable at Z\ 
and g 2 (z) = (z — z 2 ) m2 f{z) is n— times differentiable at z 2 . Define 

M-l 

g{ Z ) = f( Z ) - ^2 i C n( Z l^ Z 2)(z ~ z l) +C n {z 2 ,Z 1 )(z - Z 2 )) (z - Zi)"™" 1 (z - Z 2 
n=0 



-n-1 



where M =Max{mi,m 2 }. Then, the thesis of Proposition 2 holds for f(z) replaced by 
g(z). Moreover, if (z — zi) mi (z — z 2 ) m2 f(z) is an analytic function in Q, then the thesis 
of Theorem 1 applies to g{z). 



4. Two-point Taylor-Laurent expansions 

Theorem 3. Let SIq and Q be closed and open sets, respectively, of the complex plane, 
and Ho C O C ff. Let f(z) be an analytic function on z\ S fio and z 2 £ fl\Q, . 

Then, for z £ Q, \ f(z) admits the Taylor- Laurent expansion 

N-1 

f(z) = [d n (zi, z 2 )(z - zi) + d n (z 2 , zi)(z - z 2 )} (z - Zl ) n {z - z 2 ) n + 

7-i ( 2 °) 



^ e n (z 1 ,z 2 )(z - z 2 ) n {z - zi) n 1 + r N (z 1 ,z 2 ;z), 



71 = 
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where the coefficients d n (zi,z 2 ), d n (z 2 ,z\) and e n (zi, z 2 ) of the expansion are given by 
the Cauchy integrals 



f(w) dw 

u n {z 1 ,z 2 ) = —— -/ 

2m(z 2 - z x ) J Tl [w - x 

and 



d n (z 1 ,z 2 ) = I Vn~7""" — y^T ( 21 ) 

2m(z 2 - z\) J Tl (w - zi) n (w - z 2 ) n+i 



e n (z 1 ,z 2 ) = ^—. ^ / , {W Z ll n +1 f(w)dw. (22) 
2m J T2 (w - z 2 ) n+l 



The remainder term r (zi , z 2 ; z) is given by the Cauchy integrals 

i s _ 1 f f{w)dw N 

r N (z 1 ,z 2 ;z )=— / t 7777 7777 r (z - z\) [z - z 2 ) 

2m J Tl (w — zi)" (w — z 2 ) N (w — z) 

1 f (w- Zl ) N f(w) dw (z-z 2 ) N 
2^iJr 2 (w-z 2 ) N (w-z) (z- Zl ) N ' 



(23) 



In these integrals, the contours of integration Ti and T 2 are simple closed loops contained 
in fl\Q which encircle Q in the counterclockwise direction. Moreover, T 2 does not 
contain the points z and z 2 inside, whereas F\ encircles T 2 and the points z and z 2 (see 
Figure 6 (a)). 

The expansion (20) is convergent in the region (Figure 7) 

D ZljZ2 = {z £ Q\Q , \(z - z{){z - z 2 )\ < ri and \z - z 2 \ < r 2 \z - z x \} (24) 

where r 1 = Inf w€€ \ n {\(w - zi){w - z 2 )\} and r 2 = Inf w€ n {\(w - z 2 )(w - Zi) _1 |}. 
Proof. By Cauchy's theorem, 

, 1 f f(w)dw 1 f f(w)dw , , 

f(z) = — / ^ — / ^ , 25 

* m JF! w — z 2m J V2 w — z 

where Ti and T 2 are the contours defined above. We substitute (7)-(8) into the first 
integral above and 

1 _ Z2-Z, 1 u ^ {W- Z1 ){Z-Z 2 ) 



w — z (z — z\){w — z 2 ) 1 — u (z — z\)(w — z 2 ) 

into the second one. Now we introduce the expansion (9) of the factor (1 — u)~ l in both 
integrals in (25). After straightforward calculations we obtain (20)-(23). 

For any z verifying (24), we can take simple closed loops T± and T 2 in f2\r2 such that 
\(z-z 1 )(z-z 2 )\ < \(w — zi)(w — z 2 )\ Vu> G Ti and \ {z — Zi)(w — z 2 )\ > \(w — Zi)(z — z 2 )\ 
V w G T 2 (see Figure 6 (b)). On these contours \f{w)\ is bounded by some constant 
C: \f{w)\ < C. Introducing these bounds in (23) we see that lim^v^oo J"jv(^i, z 2 ; z) = 
and the proof follows. □ 
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(a) (b) 
Figure 6. (a) Contours T\ and T 2 in the integrals (20)-(23). (b) For z G D Zl>Z2 , we can take a 
contour Y 2 situated between and D ZlZ2 and a contour T\ in Q, which contains D Zl:Z2 inside. 
Therefore, \(z — Z\)(z — z 2 )\ < \(w — z\){w — z 2 )\ V w G Ti and \(w — Z\)(z — z%)\ < 
\(z - z 1 )(w - z 2 )\ V w G T 2 . 




(a) 4n > \zi - z 2 \ 2 , r 2 > 1 (b) 4r 1 > \z 1 - z 2 \ 2 , r 2 = 1 (c) 4r 1 > \z x - z 2 \ 2 , r 2 < 1 




(d) 4ri = |zi - z 2 \ 2 , r 2 > 1 (e) 4r x = \z x - z 2 \ 2 , r 2 = 1 (f) 4n = |zi - z 2 \ 2 , r 2 < 1 




(g) 4ri < |zi - z 2 \ 2 , r 2 > 1 (h) 4ri < |^i - z 2 | 2 , r 2 = 1 (i) 4ri < |zi - z 2 | 2 , r 2 < 1 

Figure 7. 7/?e region D ZlyZ2 defined in Theorem 3 is given by D ZljZ2 = D\ f] D 2 , where D\ is 
the Cassini oval of focus z\ and z 2 and parameter r\. On the other hand, forr 2 < 1 (r 2 > 1), D 2 
is the interior (exterior) of the circle of center Zi + (1 — r|) _1 (z 2 — Z\) = z 2 +r 2 (r 2 — l)~ 1 (zi— z 2 ) 
and radius \z\ — z 2 \r 2 /\r 2 — 1|. For r 2 = 1, D 2 is just the half plane \z — z 2 \ < \z — Zi\. The 
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shape of the Cassini annulus depends on the relative size of the parameters sjr\, ^fr^ and the focal 
distance \zi — z 2 \. 

If the only singularities of f(z) inside fio are just poles at zi, then an alternative 
formula to (21)-(22) to compute the coefficients of the above two-point Taylor-Laurent 
expansion is given in the following proposition. 

Proposition 4. Suppose that g(z) = (z — zx) m f(z) is an analytic function in $7 for 
certain m 6 N. Then, coefficients d n (zx,z 2 ) and d n (z 2 ,zx) in the expansion (20) are 
also given by the formulas: 



do{z 1 ,z 2 ) = 



(m — k — 1) 



/O2) _ 1 9 

z 2 - z x £f Q (m-k- 1)! (z 2 - zi 



(27) 



do(z 2 , Zx) 



1 9 {m) {zx) 
m\ zx — z 2 



and, for n = 1,2, 3..., 



d n {zx,z 2 )- -j— ^ U7 - ■ - — - — — x j- An+k+2 + 




(n + k)\ g (m+n-k-i)r z ^ 
k\(m + n — k — 1)! (z 2 — z\] 

f {n - k) {z 2 ) I 



k\(n - k)\ (zx -z 2 ) n+k+1 



(-1)" f ^ (n + fc-1)! g( m +"- fc )(z!) 

dnlZ2,Zlj " n! ] n 4^ fe!(m + n-A;)!(z2-«i) n+fc+1 

k fc-o ^ 

g (n + A;)! / ( "- fc - 1) (z 2 ) [ 



fc!(n - fc - 1)! (zi - z 2 )»+fc+2 



For n = 0, 1,2, coefficients e n (zx,z 2 ) are given by 



e n (zi,z 2 ) = 



(-1)' 



(n + k)\ 



9 



(m — n — k — l) 



CO 



n! 



m — n — 1 

^ fe!(m-n-A;-l)! (z 2 - zi) n+fc 

fc=0 



(30) 



Proof. We deform both, the contour Tx in equation (21) and the contour T 2 in equation 
(22) to any contour of the form Cx UC 2 contained in O, where Cx (C 2 ) is a simple closed 
loop which encircles the point z\ (z 2 ) in the counterclockwise direction and does not 
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contain the point z 2 {z\) inside (see Figure 3 (c)). Then, 



d n (zi,z 2 ) = 



dw 



2m(z 2 - zi) { J Cl (w - z 2 ) n+1 (w - z-i) n + m 
f(w) dw 
C2 {w - Zl ) n (w - z 2 ) n + 1 



Jn+m — 1 



9i{w) 



(z 2 - z x ) 1 (n + m - 1)! dw n+m - 1 (w - z 2 ) n+1 
1 d n f(w) 



+ 



W = Zi 



n\ dw n (w — zi) r 

an analog formula for d n (z 2 ,zi), and 

zi - z 2 f g(w) 



e n (zi,z 2 ) 



dw 



2vri J Cl {w - z 2 ) n + l (w - z x ) m ~ n 

1 tfn,-n-\ < w \ 

[Zl - z 2 ) 



{m - n - l)\ dw m - n ~ l (w - z 2 ) n+1 



W = Zi 



From here, equations (27)-(30) follow after straightforward computations. □ 
Remark 3. Let z be a real or complex variable and zi and z 2 (z\ / z 2 ) two real 
or complex numbers. Suppose that (z — z\) m f(z) is n— times differentiable at z\ for 
certain m £ IN and f(z) is n— times differentiable at z 2 . Define 



m — 1 



g(z) = f(z) - e n ( Zl ,z 2 )(z - z^-^z - z 2 ) n . 



71 = 



Then, the thesis of Proposition 2 holds for g(z). If moreover, {z — z\) m f(z) is an analytic 
function in CI, then the thesis of Theorem 1 applies to g(z). 
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